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Abstract 



We will introduce the biological motivation of the 7- food-limited model with variable parameters. New criteria 
are established for the existence and global stability of positive periodic solutions. To prove the existence of steady- 
state solutions, we used the upper-lower solution method where the existence of at least one positive periodic 
solution is obtained by constructing a pair of upper and lower solutions and application of the Friedreichs Theorem. 
Numerical simulations illustrate effects of periodic variation in the values of the basic biological and environmental 
parameters and how the adaptive harvesting strategies affect fishing stocks. 
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I. Time- Varying Fishery Models: Biological motivation 

Most populations experience regular or recurring fluctuations in biological and environmental factors which 
affect demographic parameters |[3|, (8), |[T0|, |[Tp|, (15], and mathematical models cannot ignore for example, 
year-to-year changes in weather, the global climate variability, habitat destruction and exploitation, the expanding 
food surplus, and other factors that affect the population growth Q, |11|, |16|, |20|- pA\ and |27|. In most 



models of population dynamics, increases in population due to birth are assumed to be time-independent, but many 
species reproduce only during a single period of the year. There are several biological parameters that can vary 
seasonally, including some cyclical changes in control parameters. For example, in temperature or polar zones 



growth frequently slows down, or even ceases in winter. Careful analysis in |21| shows that there might be a 
relationship between asymptotic recruitment and bottom temperature, i.e., stocks located in warmer waters had 
lower asymptotic recruitment. 

Consider the following autonomous model for the harvested population with size N{t) at time t : 

^^ = «("'-^ m 

where 

G(0) = r, G{K) = 0. 

We assume that G{N) is strictly decreasing, and define the intrinsic growth rate r > for w 0, the carrying 
capacity K > 0, and the effort function E. 

The linearity assumption in the logistic model is violated for nearly all populations, e.g. for a food-limited population 
Smith \25j (see 1^26 1) reported a snag in the classical logistic model, i.e. it did not fit experimental data, and suggested 
a modification of the logistic equation 



where G{N\r,K,f3) = r(l - N/K){1 + PN/K)-^, and Smith ||25J called the coefficient (1 + pN/K)'^ the 
delaying factor. 

Let x{t) = then equation Q has the form 

dx 1 — X 

— = rx Ex 

dt 1 + (3x 

To take into consideration a crowding factor , we introduce a new function 

G(x|r,/3,7) = r(l-x^)(l + /3x^)-^ (3) 

where 7 > 0. Then equation ([T]l takes a form 

dx I — x"^ ^ 

— = rx- Ex. 

dt 1 + px'y 



It is clear, if /3 = 0, the function G{x) in ([3]) is concave and G{x) — ?■ — oo as x — ?■ oo. If /? > 0, then G{x) is 
sigmoidal with — ^ < < r. 

Remark. Let Gi{x) = x{l — x'^){l + Clearly (see for example Fig 1 below) function Gi{x) has a 

unique maximum on the interval x G (0, 1), and the maximum value and the position of a critical point depend on 




Fig. 1. Different forms of function Gi{x) (Li: /3 = 0, 7 = 1; L2: /3 = 0.2, 7 = 5; L3: /3 = 4, 7 = 0.5) 



some combinations of the parameters 7 and /3. If /? = and 7 = 1 then Gi{x) is a classical symmetrical logistic 
function. 

We construct a time-varying (nonautonomous) model based on equation ([2]) by allowing r{t) and K{t) to vary, 
while maintaining constant parameters (/3,7). Similar model for unharvested population was studied in [11] under 
the assumption that K and r oscillate with small amplitudes. For an exploited marine population we introduce a 
varying effort function E = E{t). 
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No > 0. 



where 

In population dynamics the term is refereed to as the Richards' nonlinearity |26|. 

Denote 
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N{s)ds 
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where T > is the period of the system. Illustrative periodic functions K{t) and r{t) can be of the forms: 

27r(i — 



r{t) = ro 



1 + ar sin ■ 



1 + ax sin 



2lT{t 



Tk 



(5) 



(6) 



Parameters ro, Kq, a, and T can be different for r and K. We assume, however, that periods have a lowest common 
multiple, which defines the period for the system. 

For a canonical logistic equation it was proven that N decreased as the magnitude of variation in K increases, and 
that N < K irrespective to r. However, it was shown in [ 1 1 J that the effect of the environmental cycles is "a very 
model-dependent phenomenon". 

Harvesting population models with a periodic function E{t) have been studied extensively in recent years |j5J- 
1 13 1, 1 17 1- |20|. However, function E{t) is not a periodic function, in fact, it is a function of many (continuous) 
variables which fishery managers can manipulate and this function can be defined by the fishery strategies. For 



example, adaptive harvesting strategies in fisheries 1 12 1 rests on a combination of three elements: (a) deduction from 
prior knowledge of the ecosystem's components, (b) experience with similar ecosystems elsewhere, coupled with 
(c) mathematical modeling. Adaptive fishery strategies fit better with new multi-species management, with more 
emphasis on an "ecosystem approach" to sustainable fishery management. The study of the dynamical behaviors 
of the Fox harvesting models 

in periodical environments were introduced in |[5| and |[T7|. 

The paper is organized as follows. In the next section we study qualitative behavior of the solutions of a harvesting 
model in constant environments and obtain the explicit conditions for the existence of a unique positive solution 
of equation. In Section 3 for equation Q we will prove that it possess positive, bounded, asymptotically stable 
periodic solutions. In Section 4 we will investigate numerically the effects of periodic variation in the values of the 
basic population parameters r{t) and K{t) and discuss adaptive vs static harvesting strategies. 



II. Autonomous Model with Proportional Harvesting 
Consider equation (|4]) with proportional harvesting 

^ = rx[l- x^'][l + Px^]-'^ - Ex. (8) 

If /3 = then ^ is an alternative to the logistic fishing model with the Richards' nonlinearity 

dx 

— = rx(l — x'^) — Ex. (9) 



Equation ^ has a nontrivial equilibrium point 



1 - E* 



1 + I3E* ' 



XGe 

and equation ([9]) has a nonzero equiUbrium point 

XLe =</l-E*, 

where E* = E/r and < E < r. 

For the corresponding annual equilibrium harvests Yg(,{E) = Exce and Yl(,{E) = Exie, 
XGe < XLe and Yce < Yie- Notc that the maximum sustainable yield (max^y(£')) exists for E G (0,r/g). 

Let u = lnx then equations ([8]) has the following form 

du 



where 



Note 



1 - e'^'^ 
^(u) = r- . 



— = ^ (u) = - jre^ — 

du + 



and positive equilibrium 



^ 1-E* 

T+Je* 



of ([8]) is locally asymptotically stable. 



Theorem II.l. If in equation ^ we assume that r, K, /3, E and 7 are all positive constants, then 
al) If r > E then for every solution x{t) 

lim x(t) = X* 

t—>-co 

exists. 

a2) If E = then limt^oo x{t) = 1. 

Proof. The solution to equation ([8]l has an implicit form 

E* -1 + {I3E* + l)x^ = Cx"^ exp(-ra7(l - E*)t) 

with 

1 + 

C is an arbitrary constant. As t — )■ 00 

/ 1 - ii;* 
\l + E*l3 

To prove the second part of the theorem, we note that if = then the solution of ([8]) takes a form 

— l + x'^ = Cx"^^ exp(— 7art) 



thus 



lim x(t) = 1. 



III. NONAUTONOMOUS MODEL WITH SEASONAL HARVESTING 
Let US assume that all parameters in Q are continuous functions and for all t > 



I dN , ^ K^{t)-N^ 



where 7 > and /3 > 0. 



Definition IILf. We say that a positive solution N*{t) of equation (10 1 is a global attractor or globally 
asymptotically stable (GAS) if for any positive solution N{t) 

lim \N{t) - N*{t)\ = 0. 



Usually N* (t) is a positive equilibrium or a positive periodic solution of equation ( 10 1 if it exists. In general, 
we will use the following definition. 



Definition III.2. We say that equation (10 1 is GAS, if for every two positive solutions Ni[t) and N2{t) of 



equation (10 1 we have 



lim \Ni{t) - N2{t)\ = 0. 



Note that canonical logistic equation with variable parameters has been well-studied, and the questions of the 
existence and stability in this case are easily handled since the equation is solvable as a Riccati equation in a closed 
form Q- (g, flOl, ||20| and 

Theorem III.1. Assume that 

E,{t)=E{t) + ~>^ (11) 



and r{t),K{t) and E{t) are positive functions. If < A^(0) < K{Q), then every solution of equation (lOl satisfies 
the following inequaUty < N{t) < K{t) for all t > 0. 



Proof. Let v{t) = In then equation (11) has a form 

^=r{t),^ Ei{t), (12) 

where f (0) < 0. Let us prove that v{t) < for all t > 0. Suppose there exists ti > such that v{ti) = 0. Then 

dv , ^ , X 
^k=-i^i(t,)<0. 

Therefore in some interval (ti — a, ti) the function is decreasing and v{ti) = 0. But this is impossible, therefore 
v{t) < follows by N{t) < K{t). 

To prove our next theorems for the periodic models we will use p4) (see also ||5|). 

Theorem III.2. (Friedrichs Theorem). Suppose that A^) is a smooth function with period T in t for every 
N. Suppose also that there exist constants a, h with a <h such that G{t, < < G{t, a) for every t. Then there 
is a periodic solution NQ{t) of the differential equation dN/dt = G{t,N) with period T and A^(0) = c for some 
c e (a, 6). 

Theorem 111.3. Consider equation Q. For all 7 > 0,/3 > and i > we assume that E{t),r{t), K{t) are all 
positive T— periodic functions, and r{t) — Ei{t) > for all t > 0. Then there exists a positive nonconstant periodic 
solution No{t) such that 

K{fd)e^" < iVo(O) < i^(0), 

where 

1 , , r(t)-Ei(t) 
bo < - mm hm In — V — 7ttt¥t < 0. 
7 te[o,T] r(t) + /3Ei(t) 

Proof. Let v = In , then equation Q has a form 



Clearly, $(t,0) = -Ei{t) < and 



dv 1 — p^'>' 

^=r(t)- 1 ^i(t) = $(t,t;). (13) 

dt ^ ^1 + /3e^T V > y 



The last inequality is equivalent to the inequality 

1 r{ty-E^ 
7^%(t) + /3i?i(t)>'°- 

Therefore based on Theorem pil.2[ ), there exists periodic solution vo{t) of equation ( [T3] ), such that 60 < t'o(^) < 0- 
That yields the existence of the periodic solution NQ{t) of equation Q such that 

K{0)e^° < iVo(O) < i^(0). 

Theorem III.4. Suppose all conditions of Theorem ( |III.3| ) hold. Then there exists a unique periodic solution 
No{t) > of equation ^ such that < iVo(O) < A'(0) and 

lim [N{t) - No{t)] = 0. 



Proof. According to Theorem ( III. 1 1 the solution v{t) of equation (12l is negative for all t > 0. Firstly, let us 
prove that 



Assume that 



lim 'mfv{t) > —00. 



(14) 



lim v{t) 



t—^oo 



-OO. 



Then 



and 



lim inf ^ 

t-^oo at 



lim inf[r(t) - Ei{t)] > 0, 

t—>-oo 



lim v{t) = lim 



v{0) + 



dv 
ds 



OO. 



t~^oo t—>oo 

We have a contradiction. If limt_j.oo v{t) does not exist, then there exists a sequence {t„} such that v{tr, 
and 

d/i! 

= 0. 



Then equality ( [17] ) yields 



dt 



r{tn) - Ei{tn) 



-OO 



(15) 



0. 



That contradiction proves statement ( [T4| ). Suppose A'^i(t) and N2{t) are two positive solutions of equation (12) with 

< Ni{t) < K{0) i = l,2. Then 

N2{t) 

m 

are two solutions of equation (12 1. Assume that vi{t) < V2{t). Let 

v{t)=V2{t)-Vl{t). 



and 



V2{t) = In 



Then equation ( [121 ) takes the form 



dv 
'dt 



r{t) 



1 



="27 



1 



l + /3e' 



V2-f 



fi7 



Application of the Mean- Value Theorem transforms equation ( [16] ) to 

dv 



with 



where 



and 



dt 

a{t) = r{t) 



-a{t)v{t) 
1 + /3)7e'^W'^ ' 



(16) 



(17) 



Based on the last inequality 



^ [1 + /3e^(*)7]2 

/■OO 

/ a{s)ds = OO 
Jo 

Vi{t) < C{t) < V2{t). 



liminfi'(t) > — oo, 
liminf c(t) > — oo. 



Hence liminft_i.oo > 0, therefore for every solution of ( [TT] ) we have \\m.t^oov{t) = 0, or \\mt^oo[vi{t) — 
V2{t)] = 0. Similarly, if vi{t) > V2{t), then \\m.t^oo[v2{t) — vi{t)] = 0. Summing up, we conclude \mit^oo[Ni{t) — 
N2{t)] = 0. 

Remark. If Ni{t) and N2{t) are any two solutions of equation ([4]), then statement of Theorem 3.4 is true without 
the assumption that all functions r{t), K{t) and E{t) are T-periodic. 



IV. Numerical Experiments 
Fig 2 illustrates dynamics of the population for a different set of parameters a and /3. 



25 



Fig. 2. Size of the population for different models: NI: /3 = 0, 7 = 1 logistic; N2: /3 = 0.2, 7 = 5; N3: ^0 = 4, 7 = 0.5 



Fig 3 supports a well-known feature of the population models: increase in the amplitude of the carrying capacity 
K{t), defined by equation Q, yields decrease in the average size of the population, whereas, 

Fig 4 proves that the qualitative behavior of the system is unchanged fT^ by oscillations in r{t), defined by 
equation Q, alone. 



Remark. In all numerical experiments below we used parameters E{t) and K{t) that satisfy the condition (11 




Fig. 3. Change of environment. Population size for different values of oik (NI: qa' 



= 0.1, N2: an = 0.5 and N3: oik = 0.7). 



Fig. 4. Change of demography. Population size for different values of (Nl: a,. = 0.1, N2: = 0.5 and N3: = 0.9). 




Fig. 5. Change of phase shift. Population size for different values of ok (Nl: ok = 0.1, N2: uk = 0.5 and N3: ok = 0-7). 



The qualitative behavior of the system depends critically on the phase difference of the oscillations in r{t) and 
K{t). For example, if functions r{t) and K{t) have a 180° phase shift, then the relationship between the average 
population size and environmental and demographic variations (see Fig 5) is very different from the corresponding 
Fig 3. Often debated questions are the choice of the harvesting strategies and timing of harvesting. For simplicity, 
we consider static vs adaptive fishing strategies. 

Static Fishing Strategy. Consider a fishery manager who has no access to previous fishery data. He starts fishing 
(on Fig 6 curve A^l) all year with an annual quota of 12 tons and realized that in one-year period a fishstock is 
decreased significantly. Then the manager decides to shorten a fishing season, starts it in June and fishes for six 
months with a monthly quota of 2 tons (curve N2), and thereafter, he starts in September with the quota of 4 tons 
per month and fishes three months ( curve NS ), 



N1 

— N2 
N3 




5 10 15 20 25 30 35 40 45 50 



Fig. 6. Static Fishing Strategy 



Adaptive Fishing Strategy. Let a fishery manager have access (Fig 7 first graph) to the fishery data. He decides 
to fish in March because at that time the population attains its maximum. Curve N2 represents a six-month fishing, 
starting in March with 2 tons per month, whereas on curve A^3, fishing takes place with the quota of 4 tons per 
month but in a three-month period. Clearly, the latter strategy is more efficient than the static strategy (curve A^l 
on Fig 7) which represents a greater risk of depletion of the fishstock. 



Without fishing 




5 10 15 20 25 30 35 40 45 
1 

Adaptive fishing strategy 
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Fig. 7. Adaptive Fishiing Strategy 

V. Discussion 

In this paper we illustrate various effects of the varying environmental carrying capacity and intrinsic rates 
on dynamics of marine populations. Application of the food-limited model with 7-nonlinearity is consistent with 
fishery data, and supports a well-known feature of the population models: average population size decreases with 
an increase of amplitude of variation in the carrying capacity K{t), but N < K and the dynamics of the population 
must ride on this highly variable resource. On other hand, even large oscillations in r{t) alone leave the system's 
behavior practically unchanged. Note that similar results were obtained for the Fox fishery model (|7]l in |fT7|. 
However, the qualitative behavior of the system depends critically on the phase difference between the oscillations 
in r(t) and K(t). For example, if functions r{t) and K{t) have a 180° phase shift then the relationship between 
the average population size and environmental and demographic variations is very different. 
With no access to fisheries data a static fishing strategy represents a greater risk of extinction of marine populations 



under severe harvesting, compared with adaptive strategies that save fishstock in a cost effective manner. 
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